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1 Introduction 

The relativistic stellar structure equations uniquely link the equation of state of high- 
density matter to the mass-radius relation of non-rotating neutron stars. This implies 
that a measurement of the mass-radius relation or a sufficient number of individual 
points on the mass-radius relation determine the high-density equation of state. The 
upcoming gravitational-wave detectors are designed to be sensitive for gravitational 
radiation with frequencies in the range from several 10 Hz to a few kHz. Merging 
neutron stars are among the prime targets of the upcoming gravitational-wave de¬ 
tectors, since the frequencies of the orbital motion and the postmerger oscillations 
fall in mentioned frequency range (see e.g. |T] for a review). Neutron stars in bi¬ 
nary systems are doomed to merge since the orbital motion generates gravitational 
waves, which extract energy and angular momentum from the system and result in 
an inspiralling motion of the binary components. Detection rates between 0.4 and 
400 events per year are estimated for the upcoming instruments [2]. For near-by 
events with sufficient signal-to-noise ratio there are two complementary ways of in¬ 
ferring stellar properties from gravitational-wave detections. During the last cycles 
of the inspiral phase finite-size effects are imprinted on the gravitational-wave signal 
and may be extracted if the exact shape of the signal is understood (see e.g. 0). 
Considering merger simulations with a large sample of candidate equations of state, 
the outcome of a merger for the most likely binary mass range is the formation of a 
differentially rotating neutron-star merger remnant [SEIE]. Since the structure of 
this remnant is highly sensitive to the equation state, the remnant’s oscillations [7] 
are very characteristic for the equation of state HIEI- In this contribution we describe 
how the dominant oscillation frequency of the remnant can be employed to infer the 
radii of non-rotating neutron stars. We also mention the expected accuracy of binary 
mass determinations from the inspiral gravitational-wave signal, which are required 
to interpret the postmerger signal. 
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2 Radius measurements 


The dominant oscillation frequency of the merger remnant occurs as a pronounced 
peak in the kHz range of the gravitational-wave spectrum. This peak is a robust 
feature in all models which lead to the formation of a neutron-star merger remnant and 
which do not promptly collapse to a black hole. Simulations with a fixed total binary 
mass, e.g. Mtot = 2.7 Mq, reveal that the frequency /peak of the postmerger peak 
correlates very well with the radii i?NS of non-rotating neutron stars of a chosen mass 
Mns [H IS] • (The radius of a non-rotating neutron star for a given equation of state 
is a convenient way of characterizing the equation of state and thus the gravitational- 
wave emission of a model employing this equation of state.) In principle, any mass 
Mj^s can be chosen to relate with /peak- Empirically, it turns out that the 

relation between /peak and i?NS becomes particularly tight if M^s is somewhat larger 
than the mass of the inspiralling neutron stars, i.e. larger than Mtot/2 (here we 
primarily focus on equal-mass binaries). This does not come as a surprise since the 
maximum density in the remnant exceeds the one of the initial neutron stars and is 
more comparable to the densities in a non-rotating neutron star with Mjvs ~ 1.6 Mq. 
Therefore, the radius of a nonrotating NS with 1.6 Mq is particularly suitable to 
characterize the stellar structure of the massive rotating merger remnant and its 
oscillation frequencies. A more natural choice would be Mjvg = Mtot/2, which means 
characterizing the remnant by the radius of the inspiralling stars. Also in this case 
a relation between R{M^s) and /peak is found, but with a somewhat larger scatter 
compared to the use of a larger hducial neutron-star mass Mj^s [H]- 

Our procedure of inferring NS radii from NS binary mergers relies on the inversion 
of the relation between the fiduical NS radius i?NS and the dominant GW frequency 
/peak of the postmerger phase. As already pointed out in [H [5], there exist different 
relations for different total binary masses Mtot- This is illustrated in Fig. [H which 
displays simulation results for different temperature-dependent high-density equations 
of state. The temperature dependence is an important feature of these equations 
of state for the quantitative results of the merger simulations [S]; details on these 
equations of state can be found in [9]. Binaries with mass ratios different from unity 
yield to a good approximation the same /peak as the equal-mass binary of the same 
total mass. Therefore, a corresponding i?NS-/peak relation with a hxed mass ratio q 
unequal unity can be very well approximated by the relation given by the symmetric 
binary with the same total mass HI Ej- Our recipe for determining NS radii thus 
only relies on the capability of GW detectors to reveal the total binary mass from 
the inspiral phase and of determining the dominant frequency of the postmerger 
phase. The mass measurement will allow to choose the corresponding /peak — -Rns 
relation depending on M^t- Basically, this relation could be computed after M^t is 
known from a measurement by performing a number of additional simulations for the 
different EoSs with the measured total binary mass. Alternatively, the relation for a 
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Figure 1: Gravitational-wave frequency of the dominant oscillation of the postmerger 
remnant for different equations of state as function of the radius of nonrotating neu¬ 
tron stars with 1.6 Mq. Shown are simulation results for total binary masses of 2.4 Mq 
(plus signs), 2.7 Mq (circles) and 3.0 Mq (crosses). Solid lines are least-square hts to 
the empirical relations for fixed total binary masses. 


given Mtot could be obtained by a simple interpolation between two given /peak — -Rns 
relations for two different total binary masses. 

The capabilities of the upcoming GW detector network to measure the binary 
parameters of NS mergers have been explored for instance in [10]. These simulations 
have shown that for a network signal-to-noise ratio (SNR) of 20, the total binary mass 
is determined from the GW inspiral phase with an accuracy of about 1.5 per cent at 
the 95 per cent conhdence level. As demonstrated in [11] the dominant postmerger 
oscillation frequency can be measured only for merger events within about 25 Mpc. 
For such a detection, the SNR of the inspiral phase is expected to be five to ten times 
higher. Anticipating that the precision of the A/ot determination scales linearly with 
1/SNR, the total binary mass will be known practically exactly for the purpose of 
choosing the corresponging /peak — -Rns relation. 

While small asymmetries of the binary system do not change the /peak — -Rns 
relation for a given Mtot, one may want to check that the mass ratio of a detected 
merger event does not deviate too strongly from unity. For SNR=20 the individual 
masses of the binary system are recovered within ~10 per cent with an even higher 
precision for larger asymmetries [10]. In the case of a near-by merger event allowing 
a measurement of /peak, the individual component masses will thus be determined 
within a few per cent. This implies that even /peak — -Rns relations for given (Mtot, q) 
combinations can be employed to measure NS radii, which may prove useful especially 
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for systems which strongly deviate from q = 1. In summary, these considerations 
show that the binary masses of near-by NS mergers will be measured sufficiently 
accurately from the inspiral phase and will only introduce negligble errors to the 
radius determination. 

The main source of error in a neutron-star radius measurement from the post¬ 
merger gravitational-wave signal is given by the scatter in the empirical /peak — -Rns 
relation. It is not possible to predict in which way the /peak of the true equation of 
state deviates from the empirical /peak — -Rns relation. Therefore, one has to adopt 
a possible error given by the largest deviation of any theoretically possible model 
from the empirical relation. Here we implicitly assume that the sample of candi¬ 
date equations of state contains extreme models, which yield deviations from the 
fit which are larger than the one of the true equation of state. It can be seen in 
Fig. [Uthat the error of the radius measurement is of the order of about 200 meters 
for Mtot = 2.7 Mq. For larger or smaller total binary masses a determination of 
72(1.6 Mq) yields a somewhat larger error, which is understandable considering the 
discussion above about the optimal choice of M^vs. Indeed, a smaller scatter is found 
for the Mtot = 2.4 Mq data if /peak is described as a function of radii of ~ 1.35 Mq 
neutron stars |5]. Similarly, choosing 72(1.8 Mq) to describe the Mtot = 3.0 Mq data 
leads to a tighter relation than the one shown in Fig. [H Note that comparisons to the 
simulations of other groups confirm the quantitative results for the peak frequencies, 
e.g. [IS Uni HU [lU HE] (see also discussion in [5]). 

The detectability of the dominant oscillation frequency was discussed in m- 
Therein it was shown that the frequency determination in a gravitational-wave de¬ 
tection contributes only a very small error to the radius measurement. 

Finally, we note that a measurement of /peak allows a rough radius determination 
even without the precise knowledge of the total binary mass since the fits for different 
binary masses do not differ too much (see Fig. H])- With reasonable assumptions 
about the possible Mtot range, e.g. from the precisely measured chirp mass or from 
statistical arguments derived from pulsar observations, NS radii may be determined 
within km precision. 


3 Summary and outlook 

In this contribution we presented the idea to infer the radii of non-rotating neutron 
stars from the detection of the dominant oscillation frequency of the gravitational- 
wave signal of neutron-star merger remnants. We shed light on the theoretical back¬ 
ground of the underlying relations and discussed several sources of error. We point 
out that the detection of the postmerger gravitational-wave emission of two events 
with different total binary masses may be employed to estimate the emission prop¬ 
erties close to the threshold binary mass leading to the gravitational collapse of the 
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remnant. The properties at the threshold in turn allow the determination of the 
maximum mass of non-rotating neutron stars and the radius of the maximum-mass 
conhguration [9]. Also the direct observation of a prompt black-hole formation after 
a merger and the measured total binary mass may constrain the maximum mass [6]. 
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